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Q-f Abstract 

We use the vacuum expectation value of a Wegner- Wilson loop representing a fast 
,—, , moving quark-antiquark pair to derive the light cone Hamiltonian for a qq meson. 

We solve the corresponding Schrodinger equation for various trial wave functions. The 
result shows how confinement determines the meson mass and wave function for valence 
quarks on the light cone. We also parametrize the effect of the spin-dependent splitting 
for a light meson and charmonium. The correct chiral-symmetry breaking pattern for 
the pion mass is obtained due to the self-energy of the quark. 
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1 Introduction 

One of the challenges in quantum chrome-dynamics (QCD) is to solve the relativistic bound 
state problem. In the light cone Hamiltonian approach wave functions are boost invariant and 
have a well-defined probability interpretation - in contrast to the Bethe-Salpeter equation. 
But it is necessary to know the light cone Hamiltonian, to calculate reliable light cone 
wave functions. Even for the quark- ant iquark Fock space such a Hamiltonian has not been 
derived. Various approaches have been proposed to circumvent this problem. In ref. [1], 
Simula uses the usual equal-time Hamiltonian and transforms the resulting wave functions 
into the light cone form with the help of kinematical on-shell equations. In ref. [2], Simonov 
and collaborators derive a light cone Hamiltonian in a model with certain string degrees of 
freedom. More ambitious is the construction of an effective Hamiltonian including the QCD 
gauge degrees of freedom explicitly and then solving the bound state problem. For mesons, 
this approach [3,4] still needs many parameters which have to be fixed. Attempts have also 
been made to find the valence-quark wave function for mesons with a simple Hamiltonian [5]. 
A necessary input for the calculation of a two-body Fock state is a confining potential in the 
light cone Hamiltonian. For the equal-time Hamiltonian and heavy quarks the numerical 
calculation of Wegner- Wilson loops provides the form of the confining potential at large 
distances. The continuum stochastic vacuum model [6,7] allows to generalize the calculation 
of Wegner- Wilson loops from equal time to the light cone. One computes the loop expectation 
value (W[C]) in terms of gauge-invariant bilocal gluon field-strength correlators integrated 
over the minimal surface using the non-Abelian Stokes' theorem and the matrix cumulant 
expansion in the Gaussian approximation. The stochastic vacuum model is used for the non- 
perturbative low-frequency background field, and the perturbative gluon exchange is used 
for the additional high-frequency contributions. The calculation of the expectation value of 
a Wegner- Wilson loop along the imaginary-time direction gives the heavy quark- ant iquark 
potential with color-Coulomb behaviour for small and confining linear rise for large sources' 
separations [8]. 

Since the computation of the VEV for the Wegner- Wilson loop can be done completely 
analytically, also other orientations of the loop can be chosen, e.g. a loop where the quark- 
antiquark pair moves along the z-direction. By transforming to Minkowski space-time, the 
dependence of the interaction potential on longitudinal and transverse separations of the pair 
can be obtained this way. Approaching light-like trajectories of the quark- ant iquark pair, 
we have deduced in ref. [9] a light cone Hamiltonian, which contains confinement from first 
principles. In this paper, we would like to complete that work by including quark self energy 
effects from the stochastic vacuum model. To discuss chiral symmetry breaking we will 
phenomenologically include a quark wave-function renormalization and spin-spin interaction 
to make the pion mass zero for zero quark mass. We then evaluate the change of the mass 
of the pion due to a finite quark mass. Surprisingly one finds the correct chiral symmetry 
breaking behaviour in this restricted Fock space representation. 

The outline of the paper is as follows: In section 2, we review the Hamiltonian of ref. [9] 
and add quark self energy terms, which are necessary to obtain reasonable values for the 
eigenvalues. In section 3, the variational method is used to estimate the eigenvalues of the 
Hamiltonian. Section 4 is devoted to a discussion of the spin-spin interaction and wave- 
function renormalization. In that section, we derive the behaviour of the pion mass as a 



function of the current quark mass. Section 5 extends the work to heavy quark-antiquark 
systems like charmonium, where the short-range interaction becomes important. Section 6 
contains our conclusions. 

2 The light cone Hamiltonian 

We consider a Wegner Wilson loop C which has a spatial extent Rq and temporal extension 
T in four-dimensional Minkowski space-time cf. Fig. 1. It corresponds to a quark-antiquark 
pair moving with velocity (3 

f>=^ (1 ) 

cosnf^j 

where the hyperbolic angle ip defines the boost (Fig. [T]). The angle $ gives the orientation 
of the meson in the (ar, x 3 )-plane. 
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Figure 1: Configuration of the Wegner- Wilson loop in Minkowski space-time. 

For the Wegner- Wilson loop in Minkowski space-time one finds from the confining field 
strength correlators [9] 



< W[C] > 
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(2) 
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where a = 0.18 GeV 2 is the string tension given by the gluon condensate G 2 and the corre- 
lation length a. The Minkowski-geometry enters via the factor 



a 



M 



cos 2 </>sinh 2 ip. 



(4) 



This form of the loop expectation value is consistent with the analytical continuation of the 
Euclidean expectation value with a = 1 — cos 2 (j) sin 2 6 into Minkowski space by transforming 
the angle 6 — > iif>. This analytical continuation is similar to the analytical continuation used 
in high energy scattering [10-12] where the angle between two Wilson loops transforms in 
the same way. 



In addition to the contribution arising from confining field configurations there are also 
contributions from nonconfining field strength correlators which are derived in ref. [9] and 
which will be further discussed later. The exponent giving the expectation value of the 
Wilson loop can be related to the light cone Hamiltonian as follows. We define the four- 
velocity of the particles described by the tilted loop 

u, = ( 1 ,0 ± , 1 (3). (5) 

and rewrite the loop with the help of the four velocity : 

e -ig f drA^Uu, _ e ~ig f dr( 7 A - 7 /3A 3 )_ /g\ 

The line integral of the gauge potential acts as a phase factor on a Dirac wave function if) 
which splits up into a leading dynamical component ip + and a dependent component i\)_. 
For very fast quarks the mass term and transverse momenta are negligible compared with 
the energy and longitudinal momentum. In this eikonal approximation the Dirac equation 
of the leading component decouples from the small component: 

id^ + = P- t (A-)^ + (7) 

= gA~4> + . (8) 



With /3 ~ 1 the phase factor in the tilted Wilson loop integrates A and leads to a VEV for 
the loop containing P~ ot = ^(P° - P 3 )\ P ot- 

< W r [C] >= e -7(P°-P 3 )UT (g) 

Using eq. 1 one finds for the light cone potential energy arising from the confining part of 
the correlation function a term of order O(-) = 0(-p+), where P + = -jk(P° + -P 3 ) is the 
light cone momentum. 

P~ t = i= (aR Vcos(0)2 + sin(0)V 7 2 ) ■ (10) 

Terms involving transverse momenta and masses of the same order 0{-^+) are not included in 
the loop as it has been calculated. Two of these terms give the standard kinetic energy term 
of free particles, which contributes to the total light cone energy. Terms with spin cannot 
be obtained from this simplified derivation. We introduce the relative + momentum k + and 
transverse momentum k± for the quarks with mass \x. By adding the above "potential" 
term to the kinetic term of relative motion of the two particles we complete our approximate 
derivation of the light cone energy P~ 



^ + hl)P+ + 



-U^rl + xi/7 2 . (11) 



2(1/4P+ 2 - k+ 2 ) y/2 

To complete the light cone Hamiltonian, we multiply P~ with the plus component of the 
light cone momentum P + = -^(-P + -P 3 ) and use that P + jM = y2^M to eliminate the 
boost variable from the Hamiltonian. Further, we follow the notation of reference [13] and 
introduce the fraction £ = k + /P + with |£| < 1/2 and its conjugate the scaled longitudinal 



space coordinate \[2p = P + x 3 as dynamical variables. For our configuration the relative 
time of the quark and antiquark is zero. Note that M 2 = 2P + P~ . 

Hf c {p 2 ) = M 2 = ^|) + 2a^p 2 + M 2 x\. (12) 

We have obtained the light cone Hamiltonian M 2 from the confining interaction in a Lorentz 
invariant manner, because the variables £, p, k± and x± are invariant under boosts. The 
valence quark light cone Hamiltonian has a simple confining potential. The magnitude of 
the confining potential is set by the string tension a. The effective "distance " of the quarks 
is given by scale free light cone longitudinal distance and the transverse distance multiplied 
by the bound state mass. The dynamical variables are the light cone momentum fraction 

i = k + /P + (13) 

with |£| < 1/2 and its conjugate variable, namely the scaled longitudinal space coordinate 

V2p = P + x 3 . (14) 

The effective "distance" between the quarks is given by the scale-free light cone longitudinal 

distance p and the transverse distance x± multiplied by the bound state mass. Note that 

the transverse confinement scale is related to the self-consistent mass of the bound state M 2 . 

In the following we will always look for solutions which obey this self consistency condition 

without explicitly refering to it. 

The transverse momentum and the longitudinal space coordinate are represented by the 

operators 

kx = -V ± (15) 

and 

The above equation f[T2l agrees in the limit of the one-dimensional motion with the equation 
for the yo-yo string derived in ref. [13]. If only the transverse motion (p = 0) is present, then 
confinement has the usual form, which is seen by setting M « 2p. In general, the squared 
mass p 2 entering the Hamiltonian of eq. (fT2"j) contains the squared current quark mass p 2 
and a quark self energy A. 

p 2 = p 2 c + A (17) 

This self energy differentiates the constituent quark from the current quark. Light fully 
relativistic quarks surround themselves with a cloud of gluons and quark-antiquark pairs. 
We will discuss the gluon part of the self energy below. For heavy quarks we will see that 
the self energy is negligible. One should also remark that in the light cone prescription 
one can not fiddle with the zero point energy of the Hamiltonian, since the Hamiltonian 
equals the operator for the squared mass. In the conventional nonrelativistic constituent 
quark Hamiltonian one always needs a sizeable negative constant oc 2^/a to get to realistic 
values for the groundstate energy of the meson. This negative constant will reappear as a 
dynamical self energy effect in the light cone Hamiltonian. One sees that a simple kinematical 



transformation of the nonrelativistic wave functions into light cone wave functions can never 
take into account such effects, since the two Hamiltonians differ in their dynamical properties. 
For light quarks the stringy confining potential is the most important part and one can 
leave out the perturbative gluon exchange. This is not the case for the heavy quark mesons. 
The other non-confining potentials from the Abelian-like part of the correlator and the 
perturbative-gluon exchange have also been worked out from the correlation function, and 
one gets for the complete valence Hamiltonian [9]: 



#29V) = hiW) + h~ 



Yukawa 



rriQT 
M 



(a 2 + k 2 ) , 2g 2 M 2 e . . 

17TrF + 2g '- 4/3 4xr ' 



with the dimensionless variable 



p 2 + M 2 x 2 ± . (20) 



The Yukawa-potential also comes from the stochastic vacuum model, which damps the short- 
range interaction at distances r > — ~ 1/(0.77 GeV), where the long-distance confining 

2 

physics sets in. Its coupling constant |- = 0.81. The same parameters have been used 
to calculate high-energy hadronic scattering in ref. [14]. For simplicity, we have fixed in 
eq. (19) the relative weights of non-perturbative non-Abelian and Abelian-like contributions 
by k = 1.0 compared with k = 0.7 in ref. [14], i.e. we have neglected the non-perturbative 
part of the correlation function, which is non-confining. 

Let us now discuss the quark self energy corrections, which are especially important when 
the current quark masses are small. Such a self energy calculation has been discussed in the 
literature, also based on the stochastic vacuum model in two calculations [15, 16] cf. fig. 2. 
In the first version the self energy is calculated only for approximately zero current quark 
mass: [15] 

4rr 

A(/, c 2 « 0) = . (21) 

It comes out negative from the confining gluon field configurations interacting with the 
quark- field. In the second version [16], Simonov considers also heavy quarks with large n 2 
and takes into account the surrounding qq state of mass M. Then the self energy has the 
form: 

A(t) = --0(t) (22) 

TV 

with 

/■oo 

<f>(t)=t / dzz 2 K 1 (tz)e~ z , (23) 

Jo 

where 

t=(ji e + M/2)a. (24) 

As before, \x c is the current quark mass, M is the uncorrected meson mass, a = 0.302 fm 
is the correlation length of the field-strength correlator. The dependence of the self energy 
correction A (t) on t is shown in Figure [3j 

6 
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Figure 2: The left-hand diagram is a sketch of the self energy correction corresponding 
to [15]. The right-hand diagram shows the improved self energy correction corresponding 
to [16]. Crosses symbolize the color field-strength correlation functions. 
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Figure 3: self energy correction A (//;?, M) for vanishing current quark mass /i c 
function of the uncorrected meson mass M. 
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The second self energy correction of eq. (1221) agrees with the constant self energy correction 
of eq. (12T1) for fi c = and M = 0. The self energy correction is negative for light flavours and 
vanishes for heavy quarks, i.e. for heavy-meson masses M. Such a functional behaviour looks 
rather reasonable. Elimination of higher gg-gluon states produces an attractive interaction. 
In section 4 we will adress the question of chiral symmetry breaking based on the outlined 
constituent quark picture. In order to do that we will parametrize the spin-spin interaction of 
the quarks in the meson in a rather crude way. For the spin-spin interaction we have not done 
any derivations and must rely entirely on a parametrization which allows us to demonstrate 
the effect of chiral symmetry breaking in the two body Fock space Hamiltonian. So we 
cannot show how the pion becomes a Goldstone boson in the light cone picture. This is a 
question which the constituent quark model in the nonrelativistic version has not been able 
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to answer either. 

3 Variational solution of the light cone Hamiltonian 

As a first step, we evaluate the Hamiltonian for zero current quark masses p 2 = 0, with the 
quark self energy correction A(0), given by eq. (|2"T|) : 

< V\Hf c {A{0))\q >= M 2 . (25) 

We compute the vacuum expectation value of the Hamiltonian eq. (1), using a variational 
method. Simple trial wave functions factorize in a longitudinal wave function 0(£) and a 
transverse wave function ip(x±). We take the following two trial wave functions (i = 1,2), 
where the first one has the conventional form of ^-dependence: 

*«(£,£l) = <M0 ■ <p(xx) for % = 1, 2 (26) 

1 r x 2 i 

<p(x±) = -^- • ex P [ _ 2~2 ( 27 ) 



v/vrxo 
with0 1 (O = v / 6-Q-e 2 ) (28) 



1/2 



and0 2 (O = yf-Q-e 2 ) , (29) 

Xo being the meson transverse extension of the meson. 

The wave functions vanish at the kinematical boundaries (£ = ±|) which correspond to the 

limits of relative infinite longitudinal momenta in the non-relativistic description: 

% U = -\>x±) = *< U = \>*±) = 0- (30) 

A non-trivial expectation value arises in the calculation of the square root operator < 
\P | \J p 1 + M 2 x^_ | \I> >. To evaluate this matrix element, we perform a Fourier transformation 
from £-space to p-space: The wave functions have a discrete Fourier representation due to 
the finite interval in £-space £ £ [—1/2, 1/2]: 

oo 

HO = Yl &(Pn) exp {-I ■ (2fl + 1) 7T • . (31) 

n=— oo 

The conjugate variable p n take the value 

p n = (2n + 1) vr. (32) 

The normalized discrete-set coefficient functions are orthogonal 

/„ (0 = exp (i • (2n+ 1) 7T • , (33) 



fn(0r m (0^ = 5 nm . (34) 

Calculating the Fourier transform of the wave function 0i(£), for example, we have 

- /3 J!((n + 1/2) tt) 

<k(p») = fa 2^TI • (35) 

With this wave function we evaluate the expectation value of the square-root operator nu- 
merically and then approximate it by a simpler function which can be used more directly for 
the evaluation of the self-consistent mass M. The variable y = M|xj_| encodes the depen- 
dence of the confining interaction on the transverse extension of the bound state. We define 
a function g(y) 

oo 

g(y)= E 8W\/S+?. (36) 

n=— oo 

which enters the complicated expectation value 

<y p 2 + M i x 2 ± ^ = Jd 2 x ± if 2 (x ± )g(y) (37) 

and an approximation to g(y) 

G(y)=y(l-e~y/ a )+p e-y 2 / b (38) 

with 

oo 

Po= ^2 ^liP^Pn ~ 3.53. (39) 

n=— oo 

In the limit of large y, the exponentially decreasing terms in (|38|) are negligible and, because 

fif(y)=y + 0(y- 2 ). (40) 

The behaviour at small y constrains the parameters a and b 



G(l/) = Po + Z/ 2 ^--^J+0(^), (41) 

--t = ^E^^)-^ - 15505 - ( 42 ) 

a 6 2 ^ p n 

Fitting G(y) to g(y) we determine a = 6.19 and obtain from eq. (J421) 6 = 550. The 
approximate function G(y) can then be integrated over the transverse space coordinate, and 
one gets: 

2^.2 



, ,, , On \/tT M X n 



7T 



M 2 x 2 /4a 2 



Mx + 



M 3 



J r 



2a? 



Erf 



Mx c 
~2a 



- 1 



where Erf (x) is the error function. We repeat the same procedure for the other wave function, 
eq. (12"9~1) . In both cases, we get self-consistent transcendental equations for M, which can be 
solved numerically. In Fig. 2, we plot the resulting masses as a function of the transverse- 
extension parameter xq of the trial wave functions \l/j. The trial wave function \l/i leads to 
a smaller value of the meson mass, which lies in the expected range of light vector-meson 
masses. The higher mass corresponding to the trial wave function ty 2 comes about from the 
higher longitudinal momenta in this wave function. The rms-extensions v/< x\ > = Xq^ of 
the mesons can be read off from the minima of both curves. We obtain 



and 



The corresponding mass values are 



and 



xo,i = 0.8 fm 
xo,2 = 0.86 fm. 

Mi = 0.85 GeV 
Mo = 0.89 GeV. 



(44) 
(45) 

(46) 
(47) 
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Figure 4: M(x ) fc> r the trial wave functions \l/i (full line) and ^2 (dashed line). The 
Hamiltonian includes the self energy correction A (u 2 = 0) o/eq. fT2Tl ). 
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4 Chiral symmetry breaking 

The method of field strength correlators is well established and we consider the presented 
extension of this method to the light cone another sign of its significance. We would like to 
differentiate from this work the purely phenomenological spin-dependent interaction which 
we will discuss now. In general, light mesons are influenced strongly by the spin-dependent 
part of the quark-antiquark interaction. It is well-known that the form of the spin-dependent 
interaction on the light cone is not simple and the restoration of rotational invariance on the 
light cone may be nontrivial. Here, we assume that this invariance can be established and 
introduce the spin dependent interaction on the light cone by a term ^-B\Oi which does not 
depend on the transverse separation nor longitudinal coordinate. We simply do not know 
the correct dependence. It should be short-ranged in transverse space and non-local in £ 
space, otherwise the spin-splitting of mesons with non- vanishing orbital momentum, like the 
&i and bi, cannot be explained. 

In order to show the effect of chiral symmetry breaking, we must make the pion mass zero 
for zero quark mass. To have the correct p-mass, we introduce an additional wave-function 
renormalization constant Z^ on a purely phenomenological basis. Such a term may arise in a 
way similar to the Bag model [18], i.e. also on the light cone one expects Casimir corrections 
from the elimination of higher Fock states in the valence quark approximation. With these 
two parameters the theory looses any predictive power for the meson spectrum. The choice 
of the two parameters solely defines a good chiral symmetric starting point to study the 
perturbation of the chiral spectrum. It needs further investigations how the transverse and 
longitudinal kinetic energies get renormalized individually. We modify the Hamiltonian in 
the following way 

H = HIZ 1 + H s L c, (48) 



HE = ^|^ + 2^ + M^ 



m = c fw- (49) 

With Z^ = 0.4, A(t) = ~®(t) and c qq = (0.77 GeV) 2 . As before, we start with current 
quark masses equal to zero. Then the squared mass of the p-meson and Goldstone mass of 
the 7r-meson are given by: 

M 2 p =< H™ > +\c qq - (50) 

and 



4 



Ml =< H% > — A c qq . (51) 

Our goal now is to show that the valence picture on the light cone is consistent with chiral 
symmetry breaking. Chiral symmetry breaking has been a challenging aspect of light cone 
theory. It is known in equal-time theories that the vacuum is very complicated and higher 
Fock components of the quark-antiquark wave function are needed in order to reproduce 

11 



the low-energy properties of the pion correctly. An interaction of the Nambu-Jona-Lasinio 
(NJL) type leads to a quark condensate which is spread out over all space. Such condensates 
are contrary to the naive light cone picture of a trivial vacuum. The excitations of this 
condensate are massless Goldstone pions. In the light cone approach, the most developed 
calculation uses the NJL- model with a vector interaction [17] and obtains very interesting 
differences of the light cone wave function between the vector mesons and pions. In our 
framework, the complicated self energy correction A(£) of the constituent quark can give 
the correct chiral-symmetry behaviour of the pion mass. We apply the Feynman-Hellmann 
theorem [19] to the light cone Hamiltonian, which has dimension [mass] 2 

BMl ,8H, 



<9/i c \ d/j, c 

and investigate what happens to the 7r-mass squared M 2 = 0, when the current quark mass 
/i c increases to finite values fi c ^ 0. Especially one may ask whether the Gell-Mann-Oakes- 
Renner relation still holds. How can the pion mass squared vanish linearly with the quark 
mass? A naive kinetic term cannot do that because then AM 2 oc [j? c . In the Hamiltonian H 
( 149]) with A(t) we have, however, 



dM. 



d/j c 



4a I 1 \ d<f>(t) 
1/4-e/ <9M 



t=(M /2)a 



^ c =0 1 - a < x± > 



(53) 



The t-dependence of A(£) cf. eq. ( 122]) influences the expectation value of H, which is 
evaluated with ^ of eq. (|26|) . For M we take the averaged meson mass of < H^c > anc ^ 
for the transverse extension we have Xq = 0.8 fm. We get a linear dependence of the square 
of the pion mass on the quark mass with a positive slope which is related to the behaviour 
of the quark self energy when one goes from light quark systems to heavy quark systems. It 
is naturally positive, because the quark self energy starts with a negative value for the light 
quarks and becomes zero for heavy quarks. 

dM 2 

— ^ « 3.38 GeV. (54) 

When we compare this value with the one from the Gell-Mann-Oakes-Renner relation [20,21] 

M 2 //i = -2 <0| ! 9 2 |0> » 3.20 GeV, (55) 



we only find a small difference between our light cone calculation of eq. (1541) and the empirical 
value ^— of eq. (l55l) .For the absolute value of the quark condensate we took (—0.240 GeV) 3 
and F n = 0.093 GeV [22]. This result due to the self energy correction A(t) asks for further 
studies of the self energy correction in the light cone theory. Here, new possibilities are 
opening up in the AdS/QCD approach [23,24]. However, it should be mentioned that the 
quark condensate is scheme and renormalization scale dependent, i.e. eq. (155]) is very 
sensitive to small changes of the quark condensate. 
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5 Heavy quarks 



For heavy quarks, chiral symmetry is not relevant, but the short-range interaction Hy u kawa 
can be tested. In the derivation with field strength correlators, it comes about from the 
perturbative Abelian field strength correlator which dominates the short distances and fades 
out at large distances r > — « 1/(0.77 GeV). The coupling constant |- = 0.81 of this 
"Yukawa" part of the potential is taken over from the successful high energy calculations of 
hadron-hadron scattering and e-p scattering at Hera energies [14] . We can test its magnitude 
and form comparing with the mass spectrum of heavy quarks. For heavy quarks we have 



H 



QQ 
LC 



zjren. 
H LC 



+ H& + H: 



l LC 



Yukawa- 



H 



QQ 
LC 



(1/4 -e) 



+ 2a,/ p 



,« 2 i QQ - =; 



-4/3 



2g 2 M 2 e 



4nr 



(56) 



/ it en . 



m 



Hyv 



Let us take the charmonium system as an example. The charm mass is chosen as p 2 = 
1.82 GeV 2 and the self energy A(£) vanishes for heavy quarks as shown in fig. (3). The 
wave function renormalization of the heavy quark can be chosen the same as for the light 
quark, which means that the higher Fock states eliminated in the valence approximation 
have approximately the same lower threshold. In other words, the important contributions 
come from higher orbital excitations or quark gluon states which lie well above the charm 
quark mass. In Fourier space the variable r in the light cone Hamiltonian has a purely 
algebraic form without any differential operator, 



p 2 n + M 2 x 



tt2 



(57) 



Therefore, the normalized wave function ^i in the Fourier representation depends on p r 

(2n + l)n and x±. 



n=—oc 



iPn 



(58) 



%l(Pn,X±) = 4>l(p n ) -(P{X±) 



3Ji((n + i)vr) 1 
2 2n + 1 v/ttxq 



cxp 



■■ 

2x1 



(59) 



where <f>i(p n ) is given by eq. fl35l) . Because of the factorizing ansatz, we can split the 
evaluation of the "Yukawa" interaction into two steps: First, we average over the continuous 
^-dependent part of the wave function 



0"n =< ^PlHyukawal^ > z 



3n 



g l M l exp 



[tx + 27m) 2 ^g x o 



M 2 






o 



(60) 



Mx 



Erfc 



v /(7r + 27m) 2 m G x 
Mx 2 
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Secondly, we sum over the discrete Fourier components. We write the second averaging 
process in the form 



<< Hyukawa »~ 



oo 






a n 4>i{pr, 



5 

E 



On^l (Pr, 



(61) 



In practice, we truncate the infinite sum in eq. ( 16T1) by the first 11 leading terms. The 
proof of convergence of eq. ( 16TI) is given in Appendix A. As one can see in Fig. 4, the 
Yukawa interaction lowers the charmonium mass by about 220 MeV. The spin averaged 
mass of the Is — cc state comes out as < H t ^q + H Yukawa >= 3.06 GeV with the same 
wave-function renormalization factor Z^ ^0.4 for the kinetic term of the Hamiltonian as 
we used before. The size of the transverse-extension parameter for ls-charmonium is much 
smaller, xq = 0.3 fm. As one sees, the valence-quark Hamiltonian can be made also a reliable 
instrument for heavy-quark spectroscopy. Since we kept the wave function renormalization, 
the spin averaged mass for the heavy quarks is a prediction. The masses of the ip and r\ c 
mesons can be fitted with a value: 



-QQ 



/4 = 0.18GW 2 . (62) 

This value is about the same as the light quark parameter c qq -/4 = 0.15GeV 2 . The corre- 
sponding heavy light systems Dg/D and Bq/B have approximately similar splittings with 
CQq/A = 0.15GeV 2 (D) and 0.12GeV 2 (B). Thus the strength of the spin-spin interaction 
added to the square of the mass operator seems to be largely flavor blind. 
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Figure 5: M(x ) from Hlc + Hyukawa (full line) and from the Hamiltonian of eq. (IT21) 
without Coulomb interaction (dashed line). We use as a trial wave function \&i given by 
eq. (|2SD. 
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6 Conclusion 

Gluon field-strength correlators yield the most important confining interaction of the light 
cone constituent Hamiltonian derived in ref. [9]. Extending the range of applications of 
this Hamiltonian, we have calculated in this paper light and heavy meson masses and wave 
functions. Let us retrace the most important parts of our calculation. 

We have added to the Hamiltonian of ref. [9] a negative self energy correction A, which 
also has been calculated in the framework of the gluon field strength correlators [15, 16]. 
This self energy is necessary to obtain reasonable light meson masses. The improved form 
of the self energy correction varies with the uncorrected meson mass and vanishes for large 
meson masses, i.e. for large quark masses. The self energy correction comes about from non- 
perturbative binding corrections inside the meson. The dependence of A(t) on the current 
quark mass is crucial to obtain the correct symmetry pattern for the pion mass. To get 
to zero pion mass for zero current quark mass, a phenomenological spin-spin interaction is 
necessary. For charmonium the "Yukawa" part of the QQ interaction from the Abelian- 
like field-strength correlator is relevant. It is important that the purely phenomenological 
parameters do not vary much for both light and heavy systems. Our variational calculations 
are based on trial wave functions which factorize. This is not necessary. Extending the 
set of basis functions we found a small variation of the mean transverse momentum with 
the longitudinal momentum. Further work is needed to find out theoretically the spin-spin 
interaction and the wave function renormalization. 

The framework of the light cone Hamiltonian proposed here has to be seen in context with 
the successful parametrization of high-energy hadron-hadron scattering based on the same 
gluon- field strength correlation functions ref. [14]. Until now, the dipole approach uses 
only phenomenological light cone wave functions for the asymptotic hadronic states. The 
propagation of a color dipole in the nucleus demands a dynamical treatment [25]. The 
dipole propagates with an imaginary potential due its inelastic scatterings and with a real 
potential which confines the quarks on their way through the nucleus. Its Greens-function 
determines the shadowing of nuclear structure functions or heavy quark production in nuclei. 
Our approach may help to formulate this problem consistently. With this work, we converge 
towards a unified description of QCD bound states and hadronic scattering which has been 
a long term goal of QCD. 



Acknowledgments: We thank D. Antonov and J.-P. Lansberg for a reading of the manuscript 
and helpful discussions. 

7 Appendix 

A Proof of finiteness of << Hyukawa » 

To prove the convergence of the sum for << Hyukawa >> in eq. (I6TI) we will use the 
Abelian convergence test ref. [26]. This test says: Let {x n } and {y n } be two sequences of 
real numbers. If 
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^x n ^oo (63) 



71=1 



Un+1 < Vn (64) 

lim y n ^ oo (65) 

n^oo 

then ]T^° =1 XnVn converges. 

Now we apply this convergence test to << H Yu kawa >>• We identify x n = a n with a n of 
eq. P3) and y n = {<pi{n)) 2 with </>i(n) of eq. (J3SD, and split << H Yukawa » into ]T^° =1 x «^ 
and X^=i X -nV-n- 



For M^O and M 7^ oo we have < r < 1 with 

2mGvr 



r := lim {/x^, = exp 

n— »oo 



< 1. (66) 



M 
According to the Root test X]^Li x n fulfills the first condition eq. (1631 . 

Since (pi(n + 1) < 4>i(n) {y n } is decreasing and fulfills eq. (1641) . lim n ^ 00 0i(n) = implies 
limn-.oo y n 7^ 00 (eq. [65]). Therefore, conditions eq. fl63 | [64 1 1651) are fulfilled and Yl™=-oo x nUn 
converges. 
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